Spin and charge transport in U-shaped one-dimensional channels with spin-orbit couplings 
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A general form of the Hamiltonian for electrons confined to a curved one-dimensional (ID) channel with 
spin-orbit coupling (SOC) linear in momentum is rederived and is applied to a U-shaped channel. Discretizing 
the derived continuous ID Hamiltonian to a tight-binding version, the Landauer-Keldysh formalism (LKF) for 
nonequilibrium transport can be applied. Spin transport through the U-channel based on the LKF is compared 
with previous quantum mechanical approaches. The role of a curvature-induced geometric potential which was 
previously neglected in the literature of the ring issue is also revisited. Transport regimes between nonadiabatic, 
corresponding to weak SOC or sharp turn, and adiabatic, corresponding to strong SOC or smooth turn, is dis- 
cussed. Based on the LKF, interesting charge and spin transport properties are further revealed. For the charge 
transport, the interplay between the Rashba and the linear Dresselhaus (001) SOCs leads to an additional mod- 
ulation to the local charge density in the half-ring part of the U-channel, which is shown to originate from the 
angle-dependent spin-orbit potential. For the spin transport, theoretically predicted eigenstates of the Rashba 
rings, Dresselhaus rings, and the persistent spin-helix state are numerically tested by the present quantum trans- 
port calculation. 

PACS numbers: 72.25.-b,73.63.Nm,71.70.Ej 



I. INTRODUCTION 

Recent progress in the experimental techniques fabricat- 
ing semiconductor nanostructures 1 has made low-dimensional 
electronic transport one of the enduring focuses in condensed- 
matter physics. For one-dimensional (ID) systems, quantum 
wires (QWs) can be realized by growing nanowires such as 
semiconductor-based nanowhiskers or carbon nanotubes. In 
layered semiconductors, formation of QWs by confining the 
electron gas to a quasi- ID region is also possible in various 
ways, such as V-groove quantum wells, cleaved-edge over- 
growth, or atomic force microscopy (AFM) lithography. 2 The 
latter provides an even more flexible way of designing the 
shape of the confinement, and a quantum ring (QR) is one 
of the important examples. 

ID transport in QWs was previously focused on the charge 
properties. 3 A subsequent intensive investigation on spin- 
dependent transport was triggered ever since the proposal of 
the Datta-Das transistor, 4 whose underlying mechanism is 
based on the Rashba spin-orbit coupling (SOC) due to struc- 
tural inversion asymmetry. ? On the other hand, QRs provide 
a natural platform to study the Aharonov-Bohm effect 5 in 
solids. The idea of a "textured" magnetic field applied on the 
QR 6 opened the study of the Berry phase 7 in rings, in which 
the adiabatic transport plays a key role. The Berry phase ac- 
quired by the electron spin in rings was later discussed, 8 and 
investigation of the Rashba effect in QRs was subsequently 
initiated, 9 although the employed Hamiltonian at that time 
was "incorrect." After the "correct" ring version of the Rashba 
Hamiltonian was derived by Meijer et al. almost a decade 
later, 10 a series of theoretical discussion over the Rashba ring 
issue continued until recently. 1117 

So far we have been reviewing planar ID systems where 
the curvature either vanishes (QWs) or globally exists (QRs), 
whereas a general ID system may include a position depen- 



dent curvature. Quantum mechanical particle motion confined 
to a surface was first discussed by Jensen and Koppe 18 and da 
Costa, 19 regardless of spin, and was later generalized to in- 
clude the SOC effect. 20 When further restricted to a curved 
planar ID wire, da Costa proposed a linear potential term due 
to curvature, 19 which was later termed the geometric potential 
and has recently been observed in photonic crystals. 21 

Spin transport in a curved ID wire in the presence of SOC 
was recently discussed. 22-24 In Ref. 22, however, only the 
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FIG. 1 . (a) Schematic sketch of the tight-binding model for the U- 
channel. Left and right arms are labeled as regions I and III, respec- 
tively, both with N w sites. Region II is the half-ring part with N r 
sites, (b) In quantum mechanical approaches, the electron spin prop- 
agates through the U-channel via either the translation operator or 
the spin propagator, which is virtually identical to dragging the spins 
by hand along a U-shaped path in two dimensions, (c) Schematic 
sketch of the U-channel realized by AFM lithography. 
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Rashba SOC was considered; in Ref. 23, both Rashba SOC 
and the Dresselhaus (001) linear term (arising from the bulk 
inversion asymmetry of the underlying crystal 25 ) were taken 
into account, but both Refs. 22 and 23 did not consider the ge- 
ometric potential. In Ref. 24, the geometric potential was con- 
sidered but the SOC included only the Rashba term. Hence a 
more complete study of the spin transport in curved ID wires, 
taking into account both Rashba and Dresselhaus terms, as 
well as the curvature-induced geometric potential, is essential. 

Regardless of the geometric potential, spin precession due 
to SOCs along arbitrary paths was previously studied quan- 
tum mechanically, either through the conventional translation 
operator 23 or a spin propagator obtained by a properly defined 
spin-orbit gauge. 26 Despite the fact that the electron spin is in- 
deed forced to evolve through a ID path, these spin precession 
studies 23 ' 26 are built on a two-dimensional nature-there is no 
confinement. Thus how well this simple quantum mechanical 
picture can survive when a more realistic situation is consid- 
ered, such as a lead-conductor system subject to electric bias, 
has for years been a question we would like to answer. 

In this paper, spin precession patterns along a curved 
ID wire based on the previous formalisms, namely, quan- 
tum mechanical space translation [and its approximating 
result spin vector formula (SVF)] 23 and spin -orbit gauge 
method 26 will be qualitatively and quantitatively compared 
with those obtained by the more sophisticated nonequilibrium 
Green's function formalism 27 [or in ballistic systems free of 
particle-particle interaction, the Landauer-Keldysh formalism 
(LKF) 28,29 ]. Meanwhile, we will reinvestigate the influence of 
the geometric potential in curved ID transport. These are re- 
garded as our first goal. Whereas the SOCs in general depend 
on the momentum, electron spin traversing a curved ID wire 
encounters a varying effective magnetic field. This resembles 
the textured magnetic field 6 and is therefore closely related to 
the issue of adiabatic transport, which is our second goal in 
the present paper. 

For these purposes we consider a U-shaped ID channel, 
composed of two straight QWs and a half-QR in between, and 
theoretically inject electron spin from the source end and an- 
alyze the spin orientation along the U-channel down to the 
drain end. For computational concern, the U-channel is de- 
scretized into a finite number of lattice grid points, as sketched 
in Fig. 1(a). We label the left and right QWs of the U-channel 
as regions I and EH, respectively, each containing N w sites, 
and the half-ring as region II, containing N r sites. In addi- 
tion to the listed two goals, further investigation of the charge 
and spin transport properties based on the LKF will be the 
last goal. For the charge transport, the interplay between the 
Rashba and the linear Dresselhaus (001) SOCs leads to an ad- 
ditional modulation to the local charge density in the half -ring 
part of the U-channel and will be shown to originate from the 
emergence of the angle-dependent spin-orbit potential. For 
the spin transport, theoretically predicted eigenstates of the 
Rashba rings, 11-13,15,16 Dresselhaus rings, 16 and the persis- 
tent spin-helix state 30-32 are numerically tested by the present 
quantum transport calculation. 

This paper is organized as follows. In Sec. II we introduce 
the Hamiltonians and briefly the formalisms to be used in the 



transport calculations, which are reported next in Sec. III. Nu- 
merical results carrying out the above-listed three goals are 
reported, respectively, in Sees. Ill A, IIIB, and IIIC. Experi- 
mental aspects regarding the fabrication of the U-channel are 
given in Sec. IV. We conclude in Sec. V. 



II. THEORY 

In this section we will introduce the Hamiltonians to be 
used in the LKF, and review the different theoretical ap- 
proaches for the spin transport calculation. 



A. Hamiltonians 

In the following we first review and rederive the general 
form of the Hamiltonian for a continuous curved ID system, 
and then apply it to the ID ring case, which is the nontriv- 
ial part of our U-channel. We then write its corresponding 
tight-binding version of the Hamiltonian, to be used in the 
LKF calculation. Throughout we will not explicitly discuss 
the Hamiltonian for the straight parts of the U-channel since 
they are relatively trivial and well known. 



1. Continuous curved ID systems: General form 

Consider the motion of electrons confined in a ID planar 
curvilinear wire. Electrons originally in a two-dimensional 
plane are confined to a quasi- ID channel. We will derive the 
ID effective Hamiltonian in the presence of SOCs. Under the 
effective-mass approximation in solids, the Hamiltonian for 
an electron in our model is 

n = ^+ E E s ijP ^+v(v), a) 

i—x,y j—x,y,z 

where p = (j> x ,p v ) is the momentum operator in two dimen- 
sions, m is the effective mass, and cr^'s with j = x, y, z are the 
Pauli matrices. The second term is the general form of SOC 
in the Cartesian coordinate, where Sij is determined by SOCs 
linear in momentum such as Rashba or Dresselhaus (001) 
terms. Here V(r) represents the potential confining electrons 
to the quasi-lD channel. 

To obtain the effective Hamiltonian, we take the same ap- 
proach as Refs. 19 and 24. Let a(g : ) be the parametric equa- 
tion of a planar curve where q 1 is the arc length of the curve. 
The position of an electron in the plane can be written as 

r( q \q 2 ) = a ( q 1 )+q' 2 n(q 1 ), 

where ^(q 1 ) is the unit normal vector of a(g 1 ). V(r) is of the 
form: 

loo, else 
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We are going to obtain the effective ID Hamiltonian only de- 
pending on q 1 . The steps are (i) to write the Hamiltonian in 
the curvilinear coordinates q 1 and q 2 , (ii) to write an adequate 
transform of wave functions and (iii) to take q 2 — > 0. After 
these steps, the Hamiltonian will be separated into two inde- 
pendent parts regarding q 1 and q 2 , respectively. Let dl be an 
infinitesimal distance. We have 



dl 1 = gijdq l dq 3 



where gij, the metric tensor, is defined by the inner product, 
(dr/dq 1 ) ■ (dr/dqi ). We can use the coordinates q's and the 
metric tensor g % i to express the Laplacian: 



^2 , 1 d 

V Ip = r 

^fgdq 1 



(2) 



where g is the determinant of g^. Once we have SypV 3 in 
the Cartesian coordinates, we can obtain the expression in the 
coordinates qs via the transform's laws of tensors: 



(3) 



i= dq» 
1 dqi 

a" dq '\j 

, dq 1 dq 3 

^ dq'»dq' v ij ' 



where the primed symbols denote those in the new coordi- 
nates. Using Eqs. (2) and (3), we can write the Hamiltonian 
of Eq. (1) in the coordinates q's: 



H 



h 2 1 d 
2m ^fg dq 



ih^S^g 11 -^ 



3 = 1 



2m 



d 2 d d 



dq 22 dq 2 



- lh Y,S 2 yg 22 ^ + V(q 2 ), 

3=1 



(4) 



curvature of a(q 1 ). After the transform, we obtain 

n \ h 2 l d ( u d x V 



h 2 I d 2 X 



2m [dq 22 4/ 2 

3 

-ih Sijij 3 g 2 



Of 
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Taking q 2 — > except that in V(q 2 ), Eq. (5) becomes 

h 2 ( d 2 x k 2 x' 



2m \dq 2 ^ 



dq 2 



Renaming q 1 as s and deleting the terms dependent on q 2 in 
Eq. (6), we obtain the ID effective Hamiltonian 



h 2 d 2 



n 1D = -— 



h 2 n 2 



2m ds 2 8m 

-i7iS||(s) • a-^- - ^S_l(s) • o, (7) 
os 2 

where s denotes the arc length of the wire, and Sii and Sj_ are 
defined by 



S v <J = J2S lj (q 1 ,0)^, 

3=1 
3 

S ± - ( ? = ^^( (Z 1 ,0)^. 



(8) 
(9) 



i=i 



The second term in Eq. (7) is the curvature-induced geomet- 
ric potential, which was first introduced by da Costa 19 and 
was previously neglected in the literature of mesoscopic ring 



transport. 



1 We will later come back to investigate the role 



played by this geometric potential term. 



2. ID arc with SOC: Continuous form 



where for brevity the prime is neglected. We have used g l2 = 
g 2% = &2i in Eq. (4). The first two terms and latter two of 
Eq. (4) are not independent since g is a function of q 1 and 
q 2 . Given an eigenfunction ^(q 1 , q 2 ) of T-L, we have T-L^f = 
. Following Ref. 19, we make the transform xil 1 1 Q 2 ) = 
/ 1 / 2 ^ r (g 1 , q 2 ) with / = y/g = (1 — «g 2 ), where k is the 



Below we consider the Rashba SOC in an arc. In a two- 
dimensional electron gas (2DEG), the intensively discussed 
Rashba SOC ? reads 



4 -^(pv-pv), (10) 
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where a is the Rashba coupling parameter. In this case, Sij 
are 

S X x = 0, 



s„ 



!)■'■ 



a 



(ii) 



Syy — 0. 



The parametric equation for an arc can be written as R = rp, 
where r is the radius of the ring. The position of an electron is 
written as R e = (r—q 2 )p. Here, q 1 is r<f>. Using the transform 
law of tensors, we obtain S'^ in the q l coordinates from Eqs. 
(11), and thus 



(12) 



Using Eqs. (7), (8), (9), and (12), we obtain the confined Hr 
in the polar coordinates, 



S n 


= o, 


Svi 


—a 




= ~h 


021 


a 


= Ji' 


S22 


= 0. 



H.r = — ia —a 



, d 

rd<j> + 2r J ' 



where a 1 and a 2 are defined by Eqs. (3), or in the Cartesian 
coordinates, 

Hr = —— (cos 4>a x + s\xi4>a v ) 

r dm 



ia 
2r 



(cos<pa v — sin^cr^) , 



(13) 



which is in agreement with the terms given by Meijer et al. w 
The linear Dresselhaus (001) term, in a 2DEG expressed 



as 



25.33 



uW = ^( P v-pV), 



(14) 



can be derived similarly for the arc, but can be written even 
more conveniently by replacing from Hr with a — > (3, a x — > 

a v ,a v -> g x : 

Hd = -— (cos4>a v +sin(jxT x ) ^~ 
r am 



iB 

— (cos <fi(T x — sin 4>a y ) . 
2r 



(15) 



Thus the ID Hamiltonian in an arc in the presence of Rashba 
and linear Dresselhaus (001) SOCs reads 



H 



h 2 d 2 h 2 n 2 



2m ds 2 8m 
where H so is given in Eqs. (13) and (15) 



(16) 



3. ID arc with SOC: Tight-binding form 

Previously Souma and Nikolic derived the tight-binding 
Hamiltonian for two-dimensional rings in the presence of 
Rashba SOC. 14 Following their construction, here we take the 
ID limit, add the previously absent geometric potential term 
[second term in Eq. (7) or (16)] and the linear Dresselhaus 
(001) term, to obtain 

H = {U + 2t + U B ) a Y, c ] n c n 

n 

+ (Wn+lC^Cn+l + H.C.) , (17) 

n 

with the hopping matrix 

t n< _„ + i = -i <7 + i[coscj) nn+1 {t R a x + t D o- y ) 

+ sin<£ ni „ +1 (t R cr y + t D a x )\. (18) 

Here to = ti 2 /2ma 2 , a being the lattice grid spacing, is the 
kinetic hopping parameter, cr is the 2x2 identity matrix, 
£r = a /2a, to = 13 /2a are the Rashba and Dresselhaus hop- 
ping parameters, respectively, and <f> ntn+1 = (4> n + <t> n +i) /2 
is the average azimuthal angle between site n and site n + 1 
1 > ^n' see ^ e f- l^)- In the on-site potential term in Eq. 
(17), U + 2t responsible for energy band offset corresponds 
to the atomic orbital energy in the language of empirical tight- 
binding band calculation. In general U can also take into ac- 
count other local potentials, but here for convenience we will 
put U to zero. The additional term U g is the geometric poten- 
tial and can be reexpressed in terms of to as 



2„.2 



n k 
8m 



7T 
2N~ r 



to, 



(19) 



where relations k = 1/r and N r a = irr are used. Note that 
the U g term will be later considered only in the LKF, but not 
other quantum mechanical approaches. 



B. Spin transport formalisms 

Below we briefly review a set of different formalisms to 
be used to study the charge and spin transport in the U- 
channel. We will first introduce the tight-binding-based LKF 
(Sec. II B 1), for which the U-channel is precisely described 
by Fig. 1(a). That is, a ferromagnetic lead is attached to the 
left end of the U-channel, while the right lead is made of nor- 
mal metal; a bias potential difference is applied between the 
leads so that electrical spin injections from the left lead are 
theoretically simulated. Contrary to the sophisticated LKF, 
the quantum-mechanics-based translation (Sec. II B 2), as well 
as its approximating form the SVF (Sec. II B 3), and the spin- 
orbit gauge method (Sec. II B 4) are schematically described 
by Fig. 1(b). That is, we simply assume an ideal spin injected 
at the left end of the channel, drag the spin through a U-shaped 
path using either a space translation operator or a more elegant 
spin-orbit gauge operator, and then see how the spin direction 
changes along the path. 
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1. Landauer-Keldysh formalism 

The key role in the LKF is played by the lesser Green's 
function, which requires (i) a tight-binding Hamiltonian and 
(ii) lead self-energy. For (i), the Hamiltonian matrix for the 
half-ring part has been introduced in Sec. II A 3. That for the 
arm parts (regions I and III) can be straightforwardly con- 
structed from the first-quantized Hamiltonians of Eqs. (10) 
and (14) and will not be repeated here. The size of the full 
Hamiltonian matrix [22] amounts to N x N, where TV = 
2N W + N r is the total number of sites. Each matrix element 
is a 2 x 2 matrix because we are considering spin-i systems. 
For (ii), we consider semi-infinite discrete leads and summa- 
rize the self-energy expression as follows. 

Consider a ferromagnetic semi-infinite chain with uniform 
magnetization pointing along bm = (sin#M cos cf) M , sin 8 m 
sin 4> M , cos 9m)- Extending the nonmagnetic and continuous 
case from Ref. 27 to a ferromagnetic and discrete one, we 
obtain 



Sm (25) = t 2 c g M (E) , 



9m 



(E)= ^ g R {E - at M )\a;e M ){a-e M \ 



(20) 



where t c is the coupling strength between the lead and the 
central transport channel (and will be set equal to to), tut is 
the Zeeman splitting energy, the eigenkets are 34 



cr = + ; e M ) 



-id> ° M 

e Vm cos — 



'M 



e Vm sin — — 



— cos ■ 



and the retarded surface Green's function reads 



9 K (E) = =- 



A = 



1_ [ A - - A 2 , 
2t d | A - sgn AVA 2 - 4, 

E-(V + 2t d ) 
U 



|A|<2 
|A| >2' 



where td is the kinetic hopping parameter in the lead and 
will be again set equal to to in the later computation. The 
self-energy function, Eq. (20), is the only nonvanishing ma- 
trix element of the full self-energy matrices: [£l(25)] u and 
[Er(E)] nn . For our U-channel here, we will consider for 
the left lead 2m = O.lio to inject spins while for the right lead 
tju = to let the spins outflow freely. 

With both the tight-binding Hamiltonian and lead self- 
energy matrix constructed, one can construct the space- 
resolved retarded Green's function matrix 

[G* (25)] = {E [2] [H] [E L (25)] - \S R (25)]}- 1 , 

where [I] is the 2N x 2N identity matrix, [22] is the space- 
resolved tight-binding Hamiltonian matrix for the U-channel, 



and [Sx /ij (25)] is the self-energy matrix of the left/right lead. 
The lesser Green's function matrix is then obtained via the 
kinetic equation 

[G< (E)] = [G R (E)] [£< (E)] [G A (E)] , (21) 

where \G A (25)] is the advanced Green's function matrix ob- 
tained by the Hermitian conjugate of \G R (25)] and the lesser 
self-energy matrix is given by 

[S<(25)]=- {lZ P (E~eV p )}~lZ p (E-eV p )} f } 

p=L,R 

xf (E-eV p ), 

where /o is the Fermi function and eV p is the electric poten- 
tial energy applied on lead p. In our numerical computation 
we will put cVl = +eV~o/2 and eVu = — eVo/2 for a poten- 
tial energy difference of eVb, a bias parameter that is taken as 
positive (while the electron charge e = — |e| is negative), so 
that the electrons are injected from the left lead. In addition, 
we will consider a zero-temperature limit so that the Fermi 
function becomes step-like and will strictly cut the energy in- 
tegration range [see Eqs. (22) and (23) below]. 

Desired physical quantities can then be extracted from the 
lesser Green's function, Eq. (21), through properly defined 
expressions. 29 In this paper our main interest lies in the local 
charge density, 



e(N n ) = 



2iri 



E F +eV /2 



E F -eVo/2 



dETr s [G<(E)] nn , (22) 



and the local spin density, 

H/2 r E F+eV /2 



2-ki J E 



dE Tr s {a* [G<(E)] n J, 

Ep-eVo/2 UJj 



i =x, y, z 

where Ep is the Fermi level that will be set to 0.2to above the 
band bottom, [G < (E)] nn is the nth diagonal matrix element 
of the entire [G < (25)] matrix and is a 2 x 2 matrix, and Tr s is 
the trace done with respect to spin. The subscript on the left- 
hand sides of both Eqs. (22) and (23) stand for the nth site of 
the U-channel. 



2. Quantum mechanical translation method 

In the following we briefly review an earlier work done by 
some of us, 23 a theoretical method based on quantum mechan- 
ics to analyze spin precession along an arbitrary path. 

An electron spin injected at ro is described by a state ket 
| so; ro), where So labels the spin orientation, and is later 
evolved to another state ket |s; r) at position r, through 
the translation operator T (p) = cxp [ip/h ■ (r — ro)], i.e., 
|s;r) = T(p) |so;ro). In two-dimensional boundless sys- 
tems with Rashba and linear Dresselhaus (001) SOCs [Eqs. 
(10) and (14)], the eigenstates \±\<p k ) are well known (see, for 
example, also Ref. 23) and can serve as a convenient basis to 
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expand the spin state ket; 4> k is the propagation angle of wave 
vectors k±. Hence, expanding |so; ro) in terms of |±; <p k ), we 
can proceed by using / (p) |±; </> fc ) = f (hk±) \±;<f> k ): 

|s;r) = T (p) |s ;r ) 

= e ik ±-( r - r °)^|a;^)(a;^| S o;r ) 

cr=± 

= e <SA ' e^/^aj^Xtr^fclsojro), (24) 

er=± 

with fc = (fc+ + k-)/2,A9 = AfcAr = (fc+ -_fc_)Ar, 
and Ar = |r — ro|. The global phase involving k will be 
canceled in calculating the expectation value while the phase 
difference involving Afc = — 2mQ/h 2 with £ given later in 
Eq. (29) plays a key role in spin precession. For successive 
nearest-neighbor hoppings in Fig. 1(a), we simply apply Eq. 
(24) for every step and then calculate the expectation value 
for Pauli matrices to obtain the spin direction on each site, 
(S) = (H/2){3) = (h/2)(s;r\(<j x ,<7y,<T z )\s;r), starting with 
the assumed injected spin at the first site in contact with the 
left lead. 



3, Spin vector formula 

A further approximating step done in Ref. 23 (see the Ap- 
pendix therein) was to take the continuous limit, so that each 
section approaches to infinitesimal. After successive infinites- 
imal translations from injection point ro to a certain desired 
position r, the spinor overlaps carried by the final state ket 
were approximated as 

(ci|<7o)(c2|cri) ■ ■ • ((Tj +1 \<Jj) ■ ■ ■ (ct N \(T N -i) 

~ ( cr l| cr o) < 5(T2Cri • • • S aj+irTj ■ ■ ■ S aNIJN _ 1 , (25) 

where |uo) = |so;i"o) is the input, \tjj) is the shorthand for 
\o j] (f> k ) , 4> J k being the propagation angle of the jth section, 
and N — s- oo is the number of infinitesimal straight transla- 
tions from ro to r. A closed form of the state ket generalized 
from Eq. (24) can thus be obtained. Using the generalized 
state ket one obtains the SVF, 



/ 



h 

2 



— cos 9m cos ip k sin A0 + sin 9m 

— cos 9m sin ip k sin AG + sin 9m 



cos (<p k - ip° k + <t> M ) cos 2 4^ - cos (<p k + ip° k - 4> M ) sin 2 ^ 



2 

' • I n j. \ 2 Ae ( n jl \ ■ 2 Ae 

sm {<p k -<Pk + 9m> cos — ~ sin (tp k + <p%- <p M ) sm — 



V 



cos 9m cos A6 + sin 9m cos ((p k — (f> M ) sin AG 



, (26) 



with 

<Pk = ar g[(" cos (f> k - 13 sin <j) k ) + i(a sin <f> k - (3 cos (j) k )] , 

(27) 



= ?Sr- I (ds, 



H 2 



c 



C = V a 2 + P - 2a(3 sin 2<j> k {s) 



(28) 



(29) 



where (f> k is the 



The angle ip® in Eq. (26) stands for tp k (d 
propagation direction of the input |so; ro) 

For the present U-channel, the transport direction as a func 
tion of position coordinate s can be written as 



^ (s) = 



tt/2, s 6 [0, L] 

Tr/2-Tr^, se[L,L + wr] 
-tt/2, s e [L + nr,2L 



(30) 



L being the length of each arm; s runs from to 2L + nr. 
In the following we give two concrete examples to show the 
convenience of Eq. (26), one for the pure Rashba case and the 
other for the pure Dresselhaus, both with S x spin injection: 
(Wm) = (t/2,0). 

In the presence of only the Rashba SOC, we have from Eq. 
(27) (p® = <f> k (s = 0) = 7r/2, ip k = (j) k (s) and from Eqs. 



(28) and (29) A6 = 2(t R /t ){s/a). Putting these together 
with (9m, 4>m) = ( 7r /2, 0) into Eq. (26) we have 



lfl,S x inj 




(31) 



In the presence of only the linear Dresselhaus (001) SOC, 
we have ip k = —Ti,Lp k = —<j> k (s) — tt/2, and AG = 
2(tD/to)(s/a). Equation (26) then reduces to 



t / cos AG sin <p k 

( S >ks*inj = O COS AG COS (j) k 

1 1 -sin AG 



(32) 



Despite the elegant description of these SVFs, a crucial ap- 
proximation of the spinor overlaps that has been made in Eq. 
(25) deserves a further discussion before we move on. Take 
one pair of the overlap, say, between jth and (j + l)th. Re- 
call the eigenspinors in the presence of both Rashba and linear 
Dresselhaus (OOl)SOCs, 23 



1 

71 



(33) 



is given in Eq. (27). When the two sec- 



tions point along the same direction, i.e., (f> k = 4>^ k +1 , the or- 



thogonality becomes exact: (cr j+i\(Tj 



fT j+i cr 



. , regardless 
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of the type of the SOCs in the straight ID structure. Other- 
wise, the orthogonal approximation always contains an error. 
For the pure Rashba case, the overlap using Eq. (33) up to first 
order in A</> fc reads (a j+1 \aj) = (e 1 ^ 1 +<T J+1 <7 j )/2 = 



iA^/2 ■ 



, which indicates that the major er- 



ror term accumulating upon "turning" along the curved ID 
structure is proportional to the change of the angle A<fi k = 
</)k + — 0i and is therefore still moderate. In the presence of 
only the linear Dresselhaus term, the situation is similar. In 
the presence of both SOC terms, however, the error accumu- 
lated becomes drastic, which we will show numerically later. 



4. Spin-orbit gauge method 

The spin propagator can be obtained with the help of a 
spin-orbit gauge. 26 Noting that the highest order in momen- 
tum P = per in the Hamiltonian of a 2DEG with Rashba 
and Dresselhaus SOCs (both linear in p) is quadratic, one can 
define the spin-orbit gauge, 

A so = (A x , A y ) = ^ {a<jy + /3a x , -a<r x - (3a y ) , (34) 

to express the 2DEG Hamiltonian, 



njRD _ 1 ,nj2D,<ij2D 

= J_( P _£ A so> 2 

2m V c 



with the constant background potential V& = 
(m/ft 2 ) (a 2 + /3 2 ); recall that er° is the 2 x 2 identity 
matrix. Consider now the transformation operator, 



W u (r) = exp 



. so 



he 



(36) 



with the unitary property U &0 (r)U S0 (r) = er° ensured by 
the Hermitian A so = A so ^ from definition (34) and r be- 
ing the position vector of the electron displacement. Since the 
spin-orbit-interacting Hamiltonian of Eq. (35) differs from the 
Hamiltonian of the free electron gas (with a background po- 
tential Vb), 



p2 

H bee = — - V b a°, 
Am 



(37) 



only by the gauge term (e/c)A S0 in P, the following trans- 
formation is therefore suggested: 



U SO (r) p^SO (r) t = p + ^ [ A SO . r p] 



1 fie 

2 [he' 



A so -r, I A 



so 



(38) 



with [A so ■ r, P] = iftA S0 . Due to the noncommutability 
A y °] 7^ 0, the terms containing higher orders of |r|, in 
general, do not vanish, while in the small-displacement limit 
|r| w in which one has e (eft)" 1 A so ■ r < 1, Eq. (38) 



reduces to U so (r) PU S0 (r) f 
following transformation, 



;/cA S0 , rendering the 



(r) H tree U su (r) 



n 



RD 



(39) 



between the two systems, H RD and H flee . Accordingly, when 
|r| w or e (eft)" 1 A so • r <C 1 is satisfied, the free elec- 
tron gas "H free , Eq. (37), and the spin-orbit-interacting elec- 
tron gas H RD , Eq. (35), share the same eigenenergies E^. 
Their corresponding eigenfunctions, denoted by ^£ k (r)x!j ree 
and E k ( r )Xs°' respectively, differ from each other only by a 
phase factor, the 2 x 2 matrix U so (r), namely, \I/£ k (r)Xs = 
4> Ek (r)U so (r)xg fee . Here x s i s the s P m P art °f the wave 
function. Moreover, any wave function is constructed by a 
superposition of the eigenfunctions, so for any given wave 
function i/)(r)x lree in "H free , the corresponding wave function 
in U RD is ^(r)U S0 (r) x lree - 

The correspondence, which originated from gauge trans- 
formation (39), between H tree and % RD systems, allows one 
to construct the spin propagator for H RD ; to elaborate on 
this, consider an injected electron in system "H free described 
by ^inj( r )Xinj = E k C 'k^ £ ; k (r)]x in j with the initial spin state 
Xi n j and the weighting factor Ck- Without any spin-dependent 
mechanisms, this electron remains at spin state Xi n j> while 
importing A so turns on U so (r) so that the electron wave 
function in the spin-orbit-interacting system H RD can be ex- 
pressed by the gauge transformation in the form, 



(35) U S ° (r) Vtoj (r) X mj = £ C k ^ £k (r) W so (r) Xm 



(40) 



As a result, the spin polarization of the electron in sys- 
tem H RD varies spatially according to U 5 ° (r) Xi n j> an d thus 
U so (r) can be viewed as a spin propagator. 

For general applications based on the gauge transformation, 
assume an electron moving along an arbitrarily curved trajec- 
tory denoted as path c starting from one spatial point to the 
other. Divide this path into N pieces. Label the divided pieces 
(paths) by path 1, path 2, . . ., path N, sequentially (i.e., path 
i + 1 follows path i), an d let i"i denote the position vector 
of the displacement for the ith path. One can always choose 
large-enough N to have e (eft) 1 A so ■ r j -C 1 such that Eq. 
(36) can be approximately interpreted as a propagator for each 
Ti. The spin propagator along an arbitrary path c then reads 



U>° (r) = U so (r N )U so (r^_i) • • • U M (r 2 )U so (n) , 



which can be concisely written as 



W c so (r)=Pcxp(£ / 



(41) 



(42) 



where V is the path-ordering operator that orders the operator 
U so (rj) with earlier passing path to the right of the later 
U so (r i+1 ) such that Pcxp^e/ftc) / path J+1 _ path , A so -dr] = 

U so (v l+1 )U S0 (r<). 

Obviously, if both the ith and ( t + l)th paths form a straight 
line, then one has 



{v l+l +v l )=U™ {r i+1 )U* u (rO 



(43) 
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simply because only one dimension (component) of A so will 
be used, and thus the commutators that appeared in the higher- 
order terms of Eq. (38) vanish. In other words, if the electron 
moves along a straight line, i.e., path c is not curved, we have 
Ul° (r) = U so (r), namely, Eq. (42) reduces to Eq. (36). 

To study the spin evolution through the U-channel in the 
continuous limit, one can use U so (r) for parts I and III and 
Eq. (41) for part II. To be consistent with the discrete tight- 
binding model shown in Fig. 1(a) adopted in the LKF, how- 
ever, we will successively apply U so (r) for each nearest- 
neighbor hopping. 



III. TRANSPORT ANALYSIS 

Having reviewed the theoretical formalisms, we are now in 
a position to carry out our three goals of this paper. In Sec. 
Ill A, we compare the spin precession patterns calculated by 
quantum mechanical approaches with those by the LKF, or 
nowadays generally termed quantum transport. Meanwhile, 
we will examine the role played by the curvature-induced 
geometric potential based on the LKF. We proceed in Sec. 
Ill B with a detailed discussion for adiabatic and nonadiabatic 
transport regimes and connect the present paper with previous 
ones. In Sec. IIIC we discuss the anisotropic charge trans- 
port due to the interplay between the Rashba and Dresselhaus 
SOCs, and spin precession in special cases, which is equiv- 
alent to numerically testing the eigenstates of Rashba rings, 
Dresselhaus rings, and the persistent spin-helix state. 



A. Quantum mechanical approaches vs quantum transport 

1. Weak geometric potential 

Recall the geometric potential U g expressed in terms of 
to in Eq. (19). From the tight-binding Hamiltonian of Eq. 
(17), one can see that whether U g is sensible by the electrons 
depends on its competition with the energy band width 2to. 
Here we begin with a U-channel with N r = 50, yielding 
\U g \ w 2.47 x 10~ 4 t , which is hardly competitive with 2t , 
and N w = 30. 

In Fig. 2 we report the local spin densities calculated by 
the LKF under a high bias of eVb = 0.4io, spin components 
by the S VF, translation method, and spin-orbit gauge method, 
for the pure Rashba case with tu = O.lio in column (a), pure 
Dresselhaus case with = 0.1t in column (c), and a mixed 
case with (tji, to) = (0.1, 0.03)io in column (b). At the bot- 
tom of each column, the local charge density obtained by the 
LKF with both a high bias of eVb = OAto and a low bias of 
eVb = 10~ 3 io is also reported. At the top of each column, the 
spatially imaged spin vectors are from the LKF results. Note 
that the LKF-based spin densities (S' l n ) given by Eq. (23) have 
been normalized by requiring | (S n= i) | = h/2, while the spin 
components obtained by the quantum mechanical methods are 
inherently of unit norm due to the normalized state kets. 

Clearly in Fig. 2 all the spin curves obtained by transla- 
tion and by spin-orbit gauge methods are identical to each 



other. These curves further fit with those obtained by the 
LKF all quite well, except the oscillating tails that appear in 
the LKF results. These oscillations result from the nonequi- 
librium accumulation of the electron number that cannot be 
taken into account in the quantum mechanical approaches. 
With low bias the electrons behave like waves (as shown in the 
charge density curves in the bottom panels of Fig. 2), which is 
the assumption in the quantum mechanical approaches. The 
spin density curves obtained by translation or spin-orbit gauge 
method match those obtained by the LKF with low bias per- 
fectly (not shown). 

For the curves from SVFs, we use Eq. (31) for the Rashba 
case of Fig. 2(a), Eq. (32) for the Dresselhaus case of Fig. 2(c), 
and Eq. (26) for the mixed case of Fig. 2(b). In region I of 
all three cases, the SVF curves match all the others perfectly 
since in that region the orthogonality approximation Eq. (25) 
is in fact exact. Once the electron enters region II, the error 
contained in Eq. (25) for the SVF curves starts to accumulate. 
The error is still moderate in Figs. 2(a) and 2(c), but becomes 
drastic in the mixed case of Fig. 2(b), as we have remarked 
previously in Sec. II B 3. 

Comparing the low-bias charge densities in all three cases, 
one can see that an additional modulation appears in region II 
when both SOCs are present [bottom panel of Fig. 2(b)]. This 
charge density modulation stems from the angle-dependent 
spin-orbit potential, and will be explained in detail later in 

Sec. in CI. 

Before leaving for the stronger geometric potential case, we 
give a further discussion over the Rashba channel: Fig. 2(a). 
Since S x is one of the Rashba eigenstates in region I, the in- 
jected spin perfectly retains its spin direction until the half- 
ring section is reached. Then nonvanishing (S z ) component 
is induced in the curves from the LKF, and translation or spin- 
orbit gauge method since the eigenstates of the Rashba ring 
are no longer in-plane. 11-13 The spin direction described by 
the SVF, on the other hand, remains in-plane and perpendic- 
ular to the transport, since after the orthogonal approximation 
[Eq. (25)] the coplanar normal of a continuous ID channel is 
still concluded as the eigenstate. Hence the "generalized pre- 
cessionless" transport in the curved ID Rashba channel pre- 
dicted in Ref. 23 may not work well. The S x spin entering 
the half-ring region, in fact, starts to precess about the tilted 
eigenstate of the Rashba ring with spin precession length L so 
[given later in Eq. (45)], which matches exactly the period 
shown in (S z ) of Fig. 2(a). We will come back to this tilted 
eigenstate later in Sec. Ill C 2. 



2. Strong geometric potential 

Next we consider a U-channel with N r = 10 and the same 
N w . The geometric potential for such a N r is \U g \ = 0.025to, 
which is no longer negligible for the electrons. We keep the 
same figure orientation as Fig. 2. The only information added 
in Fig. 3 is the spin components computed by the LKF with 
low bias. 

The spin curves by translation and spin-orbit gauge meth- 
ods are again identical to each other and match the LKF curves 




X 10 



Coordinate s 



Coordinate s 



Coordinate s 



FIG. 2. (Color online) Spatially resolved spin components (in units of h/2) calculated by the LKF, translation method (Tran), spin-orbit 
gauge (SOG) method, and SVF in a U-channel with (N m , N r ) = (30, 50) under a strong bias of eVb = OAto and SOC strengths of (a) 
(in, to) = (0.1, 0)to, (b) (tii, £d) = (0.1, 0.03)io, and (c) (tj?, in) = (0, 0.1)to- The local charge densities obtained by the LKF for a high 



bias of eVb = OAto and low bias of eVo = 10 
LKF results. 



to are shown in the bottom panels. Spin vectors in the top subplots are based on the high-bias 



with high bias well. SVF curves this time become rather poor 
after entering the half-ring region since there the change in the 
direction A<fi k upon every hopping is no longer small and the 
spinor overlap approximation (25) hardly applies. 

For the spin curves obtained by the LKF with low bias, the 
effect of the geometric potential U g can now be seen. The 
injected spin previously parallel to the internal magnetic field 
direction bm is here reversed due to the reflection off the U g 
potential well and a certain matching condition between the 
Fermi wavelength and the arm length L within region I. Shift- 
ing the Fermi energy Ep, changing the length L, or putting 
different SOC energies [such as Fig. 3(b)] will make the re- 
versal of the spin direction disappear. This is also why, in the 



strong -bias regime, where a larger range of contributing states 
are integrated, the reversal does not show. 

The huge difference due to the stronger U g shown in Fig. 3 
is hence only a special case: The reflection off the U g well and 
the length matching happen to make the opposite spin state fa- 
vored upon injection. Here we conclude that the role played 
by the geometric potential is merely a rather weak potential 
well that can be possibly sensed by the electrons when the 
bending of the ID structure is severe, and that even if U g is 
sensed, it serves simply as a potential well, which becomes 
crucial only in the linear transport regime with a certain par- 
ticular matching condition between the Fermi wavelength and 
the channel size. 




Apart from the spin behavior, the potential well nature of 
the U g can be clearly identified by comparing the low-bias 
charge density curves shown in the bottom panels of Figs. 2 
and 3. In the previous weak U g case, the electron density 
distributes like a standing wave from the source to the drain 
ends, while in the present strong U g case the electron wave is 
disturbed by the central geometric potential well in the half- 
ring part. 



B. Adiabatic and nonadiabatic transport regimes 

Comparing further the previous two U-channels, one can 
see from the Rashba cases that the spin initially injected at one 
of the Rashba wire eigenstates may or may not follow the local 



eigenstate throughout the U-channel; see Fig. 4(a). Certainly 
the key lies on the half -ring part, where the larger the number 
N r is, the easier the spin can follow the eigenstate, but the 
strength of the SOC is also an important factor. 12 

To make a quantitative investigation, we first define the fol- 
lowing spin-flip ratio, 

■f^-^i^l)^^ 00 ^ (44) 

where S l and Sr are the average spin direction of the left and 
right arms, respectively, computed by the LKF In the case of 
the Rashba U-channel with S x injection, the spin is flipped 
from +S X to — S x when reaching the right arm, if the local 
eigenstate is strictly followed. Thus the definition of Eq. (44) 
helps us quantify how well the local eigenstate is followed: 
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FIG. 4. (Color online) (a) Top view of the spin density vectors shown 
in Figs. 2(a) and 3(a) with high bias, (b) Spin flip ratio f r as a func- 
tion of N r with (tn, to) = (0.1, 0)io and the output spin polariza- 
tion based on Ref. 22; the nonadiabatic curve is given by Eq. (47). 



whether a change in direction or a shrink in the magnitude re- 
duces the spin-flip ratio. Note that Eq. (44) is also valid for 
pure Dresselhaus U-channels, provided that the injected spin 
has to be oriented along ±y, due to the 180° turn of the U- 
channel. That is, as long as the spin transport is sufficiently 
adiabatic, the injected spin is able to follow the local eigen- 
state so that the spin is flipped after passing through the half- 
ring region. 

Let us first fix the Rashba SOC strength as tjt = O.lio 
but change the half-ring from small radius to larger ones, as 
shown in Fig. 4(b), where a clear jump at about N r = 16 is 
observed. Within N r < 16 the spin-flip ratio is nearly zero, 
showing that the spin can hardly follow the local spin eigen- 
state when entering the half-ring region. At right side of the 
jump, f r increases to 100% and then exhibits a resonance- 
like oscillation below the maximum value, in close analogy to 
Ref. 22 and similar to some of the results reported in Ref. 24. 
The oscillation period of about 32 corresponds to a distance 
for the spin to complete a 2ir of precession angle under the 
Rashba SOC, i.e., two times the spin precession length, 28 



L so /a = 



7T to 
2*h" 



(45) 



Here we have 2L so /a = lOir. This oscillation period can be 
well described by Ref. 22, which we will discuss later. 

The jump of f r and the wavelength of the resonance-like 
oscillation depends on the SOC strength. Hence we next vary 




FIG. 5. (Color online) Spin-flip ratio f r defined in Eq. (44) as a 
function of tj? and N r , free of to. The white dashed line given by 
Qr — 1 divides the transport into nonadiabatic (left-bottom) and 
adiabatic (right- top) regimes. Horizontal dashed line corresponds to 
Fig. 4(b). 



both Rashba strength and the site number of the half-ring, and 
make the plot for f r as a function of tR and N r in Fig. 5. The 
f r pattern is clearly divided into two regimes that can be per- 
fectly described by the Qn = 1 curve, which was inspired by 
Ref. 12. The adiabatic condition was previously argued as 8,9 
Q ^s> 1, where Q = Qb + Qr includes the contribution from 
external magnetic field and the Rashba field. In our analysis 
no external magnetic field is applied, and the adiabatic condi- 
tion reads Qr 3> 1. The definition of Qr from Ref. 12 is here 
reexpressed in terms of our tight-binding parameters, 



Qr = -^N r , 

7T t 



(46) 



which implies that the increase of either N r or tR brings the 
transport regime to adiabatic. Therefore the criterion that 
Qr ^ 1 preserves the transport in the adiabatic regime is 
well agreed. Furthermore, by a comparison with Eq. (45), 
the meaning of Qr given by Eq. (46) is transparent: Qr = 
N r / (L so /a) , i.e., the number of precession half -periods that 
the spin can complete within the half -ring. Hence the con- 
dition Qr = 1 from our notation corresponds exactly to the 
arc length of the half-ring that matches one spin precession 
length L so . The mathematical criterion for adiabatic trans- 
port, within which the electron spins are able to follow the 
local field, means that the electron spins have to be able to 
complete at least an angle of ir of precession within the half- 
ring. One might attempt to extend this interpretation to an 
arbitrary expanding angle of an arc, such as the curvilinear 
QW considered in Ref. 24, but further examination is left here 
as a possible extending work. 

Finally, we compare our result with that of Ref. 22, which 
is shortly reviewed in the following. Trushin and Chud- 
novskiy previously considered also a U-shaped ID channel 
with Rashba SOC and solved the transmission problem by 
matching the boundary conditions. 22 A spin polarization was 
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defined as P = (j + — + where j ± is the prob- 

ability current of the ± eigenspin components. The polarized 
wave occupying the + Rashba eigenstate was assumed as the 
incoming state so that Pj n = 1 and the spin polarization for 
the outgoing wave P out is the main quantity of interest. If 
the injected spin remains at its local eigenstate, P ut = 1 is 
expected, which is the adiabatic limit. Oppositely, a strong 
nonadiabatic limit leads to a simplified expression, 22 



cos ttQr\ 1 + Q 



r R 2 



(47) 



which has been translated to our tight-binding language. As 
shown in Fig. 4(b), the oscillation matches our result; the 
curves in the low 7V r region (Qr small) match particularly 
well. For larger Qr, Eq. (47) approaches cos(itQr) with the 
oscillation period AQr = 2 = AN r /(L so /a), which well 
describes the period of AN r = 2L so /a in Fig. 4(b), in agree- 
ment with our previous discussion. For general Q r, P out can 
be computed following their results and is plotted also in Fig. 
4(b) (with t R = O.lto and Ep = 0.2£ )- Overall, the oscilla- 
tion behavior in both the nonadiabatic and general cases from 
Ref. 22 agrees with our result. 

In the above discussion, we have injected an S x spin, which 
is the eigenstate of the Rashba wire. The fact that the eigen- 
states in the ring differ from those in the wire by a tilt an- 
gle is the origin that a spin starting with S x in the U-channel 
can never perfectly follow the local eigenstate in the half-ring 
part. In principle, when the spin direction happens to match 
the ring eigenstate when the electron is just about to enter the 
ring, the local ring eigenstate can then be well followed. In 
the following section, we will show that these special cases 
do exist, provided that the length L and the orientation angle 
of the injected spin are precisely designed. 



C. More on charge and spin transport 
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FIG. 6. (a) Formation of the charge density modulation in a 
(N w ,N r ) = (50,100) U-channel with eVo = 10~ 3 io, t D fixed 
at 0.02to, and ta varied from to to- (b) ID Schrodinger problem 
for a ID linear chain subject to the spin-orbit potential model. The 
top, middle, and bottom panels account for the zero SOC, single type 
of SOC, and mixed type of SOCs, respectively. Parameters used here 
are identical to the U-channel shown in (a). In addition to regions I, 
II, and III that correspond to those in the U-channel, the left (L) and 
right (R) leads are also labeled, (c) Local charge density (eN) for the 
U-channel and the electron correlation function G™ oc \ip (s)| 2 for 
the linear chain, both considering tR = to = 0.02to- Calculations 
here are all with Ef = Q.2to. 



The last goal to be carried out here is to reveal some of 
the interesting transport properties, regarding both charge and 
spin, based on quantum transport calculations. 

1. Charge density modulation: Emergence of spin-orbit potential 

As previously remarked in Sec. Ill A 1 [or specifically the 
bottom panel of Fig. 2(b)], an additional modulation to the 
low-bias charge density in the half-ring region appears when 
both terms of SOCs are present. In Fig. 6(a), we show the 
formation of this charge density modulation in a (N w ,N r ) = 
(50, 100) U-channel by fixing trj = 0.02io and varying from 
tR, = to t R = t D , with eVo = 10~ 3 t and E F = 
0.2io. Clearly, the modulation appears only when £r*d 7^ 
and reaches its maximum when both SOCs are of the same 
strength. This modulation was similarly obtained in a recent 
study of the anisotropic spin transport in mesoscopic rings, 35 
but the origin there was not clear. In the following we pro- 
vide a simple quantum mechanical picture to account for this 
modulation. 



Recall the anisotropic spin splitting, Eq. (29). By solving 
for k± from Ep = h 2 k"±/2m ± (k± and defining the Fermi 
wave vector as kp = (k + + /2, we have 




where 

tea {<j> k ) = ^ = y/t 2 R +t 2 D -2t R t D sin2^ fe . (49) 

The fact that the Fermi wavelength = 2-njkp is forced 
(when t Rtp> ^ 0) to be modulated upon changing propagation 
angle <fi k can be mapped to a linear ID chain with a position- 
dependent local potential. Defining 

Eq. (48) then reads kpa = -J (Ep — V so ) /t , as if the elec- 
tron were propagating in a ID linear system subject to a poten- 
tial V so , i.e., the electron is governed by (p 2 /2m + V so ) = 
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Eip. This interpretation is exact when V s0 is a constant poten- 
tial; when V so is position dependent but weak compared with 
Ep, the argument is still a good approximation. 

We therefore consider a ID linear chain subject to, together 
with the geometric potential, the full local potential, 



V so (+tt/2) , s e [0, L] 
{V so {<P k ) + U gi sE[L,L + irr] , (51) 

^ V s0 (-tt/2) , s E [L + tit, 2L + irr] 

where cj) k as a function of s is taken identical to Eq. (30). The 
ID Schrodinger problem subject to the potential U(s) given 
by Eq. (51), [-h 2 d 2 /2mds 2 + U(s)] ip(s) = Eij)(s), can be 
analytically cumbersome due to the irregular shape of U (s), 
but can be easily solved by the quantum transport formalism 
introduced in Sec. II B 1 at an even lower cost. The lead self- 
energy, Eq. (20), with tm = is taken for both left and right 
leads (with potentials +eVb/2 and — eVo/2, respectively) to 
simulate the incoming and outgoing waves. The squared norm 
of the wave function in such an equilibrium problem corre- 
sponds to the electron correlation function G n (E) (equiva- 
lent to —iG K ) that can be obtained from f (E) A (E). 21 The 
Fermi function fo will be taken as unity, concerning the cur- 
rently assumed zero temperature, and the spectral function 
A can be obtained from A = G R TG A , where G R and G A 
are the retarded and advanced Green's function of the lin- 
ear chain, respectively; the broadening function is given by 
r = — Note that we have turned off the contri- 

bution of the right lead to the broadening function to suppress 
the inflow of the particles from the right leads. This is equiv- 
alent to setting the wave function in the outgoing region as 

^( S )L>L+,rr 

most quantum physics textbooks. 



,+iks ^ yyjjjcjj j s a usua i assumption taken in 



The potential profile U (s) together with G n oc \ip (s)| is 
plotted in Fig. 6(b) for three situations: zero SOC, single type 
of SOC, and mixed type of two SOCs. The last case clearly re- 
sembles the charge density modulation in the U-channel with 
tiitD 7^ 0, and the present spin-orbit potential model seems to 
work well. Thus the modulation of the electron density pro- 
file simply reflects the position-dependent spin-orbit potential, 
Eq. (50), that is usually small compared with Ep. Indeed, in 
Fig. 6(c) we further compare the local charge density {eN) in 
the U-channel with the electron correlation function G n in the 
linear chain. The difference between them is only up to a tiny 
phase shift. Detailed parameters used here in Fig. 6 are given 
in the caption thereof. 

When \V so \ is close to Ep, either by strengthening the 
SOC parameters or lowering the Fermi energy, the phase shift 
grows, but the G n calculated for the linear chain and the {eN) 
calculated for the U-channel still are similar in their shapes 
(not shown). The present model hence works equally well to 
explain the charge density modulation, which we conclude as 
originating from the emergence of the angle-dependent spin- 
orbit potential when tptp, ^ 0. 



(a) t R ^0,t D =0 (b) t R = 0,t D ^0 (c)t R =t D 


i 




Coordinate a Coordinate s Coordinate s 

FIG. 7. (Color online) Spin precession in a (N w , N r ) = (50, 50) 
U-channel with (a) (t R ,to) = (-K/N w ,0)t o , (b) (t R ,t D ) = 
(0,w/N w )to, and (c) t R = to = {it/V2N w ) to. The injected 
spin is oriented in (a) as the Rashba ring eigenstate Eq. (52), in (b) as 
the Dresselhaus ring eigenstate Eq. (56), and in (c) as the persistent 
spin helix eigenstate Eq. (57). Note that in the imaging of the space- 
resolved spin vectors, only half of the vectors are drawn in order for 
clarity. 



2. Spin precession in special cases 

In this subsection we discuss spin precessions in three spe- 
cial cases. In the first two cases, only one type of SOC is 
considered, and we inject the spin oriented as the theoretically 
predicted eigenstate for the ring (which is different from that 
for the wire), and adjust tRortp such that the length L equals 
exactly two times the spin precession length L so . The injected 
spin arriving at the half-ring returns exactly to the eigenspin 
direction, such that the previously derived tilted eigenstate, 
e.g., the eigenspinor for clockwise-propagating I eigenspinor 
in the Rashba ring from Ref. 12, 



xi- (0) 



sm (7fi/2) 
e**cos( 7fl /2) 



cos[(tt - j R )/2] 
e^sin[(7r- 7fl )/2] 



(52) 

with the tilt angle 7^ = tan -1 Qr, can be numerically ex- 
amined. Note that here (f> is the azimuthal angle in the polar 
coordinate (0 = defined along +x axis), rather than the 
propagation angle qb k . 

To account for the two extreme cases of pure Rashba and 
pure Dresselhaus using one single formula, we write the 
clockwise-propagating I eigenspinor as 



X 



1(0) 



cos[(7r - 7 )/2] 
e^ sin[(7r-7)/2] 



(53) 



where tp is given, similar to Eq. (27), as 

ip = arg [t r cos cj) — tr> sin 4> + i (t r sin <f> — tp, cos < 
The tilt angle in Eq. (53) is given by 



7 = tan 



Ql + Q 2 



D 



(54) 



(55) 
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where Qjj is defined similar to Eq. (46) as = 
(2/n)(tD/to)N r . Equation (53) clearly recovers the Rashba 
ring case of Eq. (52) since (p\ t _ Q = <j> from Eq. (54) and 
7| t =0 = Jr from Eq. (55), can cover the tR = Dressel- 
haus ring case, 

i _ ( cos[(^ - 1d )/2] \ 
Xd - ~ { e -^/ 2 +^ sin[(7r - 1d )/2] ) ' P<3) 

where j D = tan _1 Q_D, but does not apply for general 
tfjtu 7^ cases. 

The last case is tR = to, corresponding to the persistent 
spin-helix in 2DEG, 30 ~ 32 the main feature of which is the fixed 
eigenspin directions. In the present case, the I eigenspinor is 

Xi=n = ^{ j J, (57) 

corresponding to Ek = ti 2 k 2 /2m — (k. Whether this eigen- 
state, valid for 2DEG, still works in the curved ID system, is 
what we are about to answer. 

(a) Rashba ring eigenstate. We begin with try = and in- 
ject Xr- (0 = 7r ) as given by Eq. (52) in a U-channel with 
(N w ,N r ) = (50, 50). We tune t R /t = tt/N w w 0.063 such 
that N w = Trto/t,R = 2L so /a [see Eq. (45)] ensures the re- 
turn of the injected spin to its initial spin direction, which is 
the eigenstate of the ring at cjy = tt, after going through region 
I. As shown in Fig. 7(a), the spin entering the half-ring region 
remains perfectly in the eigenstate. (Note that we have cho- 
sen another view angle to focus on the half-ring part; injec- 
tion conditions remain the same as in previous discussions.) 
The curves of (S l ) within region II can be well described by 

(X R - i4>) \v\Xr- (4>)) [given below in Eq. (58)] with <j> run- 
ning from 7r to 0. The validity of the previously derived tilted 
eigenstate in Rashba rings is hence numerically verified. 

Note the subtle difference between the special spin injec- 
tion here and in Sec. IIIB, where we injected an inplane S x : 
an eigenstate of the wire. The precise design of the length 
L and the orientation of the injected spin allows the spin to 
stay perfectly in the tilted ring-eigenstate. The x-component 
of spin with such a precise design can always be flipped, but 
should be regarded as a special situation. 

(b) Dresselhaus ring eigenstate. We continue with Ir = 
and inject xb- = 7r ) given by Eq. (56). In this case we 
similarly have to /to = tt/N w . Again the spin arriving at the 
half -ring enters its eigenstate and remains so until leaving the 
ring, as shown in Fig. 7(b). The spin components (S l ) within 

region II can be well described by (xb- (0) l^|Xr»- (0))> an d 
the validity of the Dresselhaus ring eigenstate is also numeri- 
cally verified. For both Figs. 7(a) and 7(b), spin components 
from Eq. (53), 

(sin (tt — 7) cos ip \ 
sin (tt- 7) sin <p , (58) 
cos (tt — 7) J 

describes the (<S l ) curves in region II well. 

We remind the reader here that the eigenstate given by Eq. 
(53) is intended only for the two extreme cases of tRm\ 7^ 



0, tij^R) = discussed above, although a solution given in 
Ref. 16, similar to our Eq. (53), was claimed to be valid for 
rings in the presence of both Rashba and Dresselhaus terms. 
The simple reason why the form of Eq. (53) does not apply 
for general cases of tRtrj ^ is that the tilt angle 7 [Eq. (55)] 
does not recover tt/2 when the tR = try 7^ persistent spin 
helix state is reached, which is true as we will next numeri- 
cally show. 

(c) Persistent spin-helix eigenstate. We proceed by con- 
sidering tR = to = (tt /y/2N w )to m 0. 044to, keeping the 
size of the U-channel unchanged. The injected spin state is 
oriented as Xr=d> given in Eq. (57). As expected, the in- 
jected spin stays at this eigenstate in region I, as shown in Fig. 
7(c). Somewhat surprisingly, however, the injected spin re- 
mains precessionless throughout the whole U-channel, even 
in the half-ring region. Therefore, the persistent spin-helix 
eigenstate, originally derived for a 2DEG, is equally valid in 
straight wires and curved rings. This is in sharp contrast to 
the pure Rashba and pure Dresselhaus cases, for which the 
eigenstates in wires and in rings are different. 



IV. EXPERIMENTAL ASPECTS 

The U-shaped ID channel theoretically discussed in the 
present paper can be experimentally prepared by using AFM 
lithography, 1 i.e., local oxidation 2 written by an AFM tip on 
the sample. The oxide lines turn out to completely deplete 
the 2DEG underneath, and can hence confine the electron gas 
in the desired nanostructure. A schematic sketch of the U- 
channel based on this technique is suggested in Fig. 1(c). Ac- 
cording to the present fabrication ability (see, for example, 
Ref. 36), however, the ring radius is of the order of 100 nm, 
and the induced geometric potential is rather weak: 10~ 2 
meV for GaAs-based quantum wells. To focus on the effect 
of the geometric potential, a single turn as in our U-channel 
is not enough; a series of geometric potential wells such as a 
sinelike waveguide similar to the design reported in a recent 
experiment on photonic crystal 21 may give rise to a resonance 
that could potentially be measured. 

The spin injection assumed here may be realized either 
electrically or optically. The former requires a ferromagnetic 
source contact and may further complicate the sample fabri- 
cation and even the transport properties. The latter, optical 
spin injection, has been mature in generating spin packets that 
can be electrically manipulated. 37 Regarding the U-channel 
sketched in Fig. 1(c), the adiabatic-nonadiabatic spin trans- 
port discussed here may be experimentally tested by optically 
pumping at the source a spin packet that can be electrically 
dragged to the drain end by applying a bias voltage between 
the source and the drain contacts. Optical spin detection of the 
spin packet at the drain end shows whether the spin is reversed 
(adiabatic), is decayed (spin relaxed), or remains (nonadia- 
batic), compared with the injected spin direction. The laser 
spot size for the optical spin injection or detection, typically 
of a few hundreds of microns, may impose a corresponding 
limit on the design but could be possibly overcome by hard 
masks. A top gate covering the U-channel may control the 
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Rashba SOC strength and switch the transport regimes 
between adiabatic and nonadiabatic, provided that the effect 
of the Dresselhaus term is well treated. 

Experimental proof of the interesting charge and spin trans- 
port properties discussed in Sec. IIIC are also expected. In 
particular, the charge density modulation in the presence of 
both Rashba and linear Dresselhaus (001) SOCs discussed in 
Sec. Ill C 1 requires measurement on the local charge densi- 
ties only and should be possible. The profile of the charge 
density modulation simply reflects the angle-dependent spin- 
orbit potential and hence determines the type of the SOCs: flat 
for the single type of SOC and sinelike for the mixed type of 
SOC. Note that in our discussion of the U-channel, the nature 
of sin 20 dependence of the spin-orbit potential provides the 
modulation in the half-ring with one period of the sine func- 
tion. In the case of a full ring, we expect two periods then. 

An alternative to preparing a ID channel is the V-groove 
QW based on electron beam lithography, 1 but so far applica- 
tion of this technique to curved ID QWs has not been seen. 

V. CONCLUSION 

In conclusion, we have rederived the Hamiltonian for a 
curved ID structure in the presence of SOC. Applied to the 
ID ring, the Hamiltonian is not only consistent with the pre- 
viously proposed proper Hamiltonian for a Rashba ring, 10 but 
also contains a curvature-induced geometric potential, which 
was first derived in Ref. 19, but less discussed in the litera- 
ture of ring issues. The U-shaped ID channel is further taken 
to be a specific example to investigate the role of this geo- 
metric potential, as well as to compare the spin densities ob- 
tained by the LKF with the previous quantum mechanical ap- 
proaches. Both translation 23 and spin-orbit gauge 26 methods 
mostly agree with the LKF results, even though the underlying 
assumption is rather simplified: to assume an ideally injected 
spin at the initial position, and the technique is rather artifi- 
cial: to drag the injected spin by operating quantum operators. 
Whether the SVF, 23 a further approximating result from the 
translation method, may work well or not depends on whether 
the orthogonality approximation [Eq. (25)] is valid enough: 
exact for a straight ID structure, of moderate error for a ring 
with a single type of SOC, and poor for a ring with a mixed 



type of SOCs. 

The influence of the geometric potential taken into account 
in the LKF calculation is shown to be sensible only when the 
turn of the channel is sharp and the transport is under low bias. 
Overall the role played by the geometric potential is moderate, 
just like a local potential well, and can be drastic [such as the 
reversal of the injected spin state shown in Fig. 3(a) or 3(c)] 
only when a certain resonance condition is reached. 

We have also discussed the spin transport in adiabatic and 
nonadiabatic regimes. In addition to the increase of the ge- 
ometric potential when making the turn sharper by reducing 
the number of site N r in the turning part, the transport be- 
comes nonadiabatic since the change of the local eigenstate 
becomes rapid. The spin transport shows adiabatic behavior 
when the turn is smooth and the SOC is strong enough, which 
agrees with the previously stated adiabatic condition. 8 ' 9 ' 12 We 
have also compared our results with a recent similar work by 
Trushin and Chudnovskiy, 22 which showed good agreement. 

The last part of the numerical results revealed interesting 
charge and spin transport properties. For charge transport, 
the interplay between the Rashba and linear Dresselhaus (00 1 ) 
SOCs leads to anisotropic spin splitting, and hence an angle- 
dependent Fermi wavelength. Charge transport in a curved ID 
chain subject to a modulated Fermi wavelength is therefore 
mapped to transport in a ID linear chain subject to a position- 
dependent potential. We have shown that the charge density 
modulation that appears only when tut u 7^ in the half-ring 
part of the U-channel can be well explained by the spin-orbit 
potential model. For spin transport, we have shown spin pre- 
cession patterns in three special cases, which are equivalent 
to numerically testing the validity of the previously predicted 
tilted eigenstates of the Rashba rings and Dresselhaus rings, 
as well as that of the persistent spin-helix state. 
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